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ABSTRACT 


The purpose of this thesis is to model digital communications systems in the time 
domain using MATLAB Simulink and the Communications Toolbox as well as to 
determine and verify system performance in the presence of additive noise and co- 
channel interference. While the theoretical results are available for the effect of wideband 
gaussian noise on the performance of digital communications systems, determining the 
performance of a system in the presence of noise = co-channel interference is best done 
by computer simulation. Time domain modeling allows the visualization of the 
communication signal at various stages, and “Monte Carlo” type simulations establish the 
bit error rates under realistic conditions of noise and co-channel interference as well as 
for different transmitter/receiver parameters, different channel parameters, and different 
types of interfering signals. It should be noted that the simulation does not account for 


system non-linearities. 
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I INTRODUCTION 


A military site or a platform typically has a large number of communication 
systems operating simultaneously in both transmit and receive modes. The interference 
of one system with another depends on available frequency allocations, the physical 
separation of transmit and receive antennas, non-linearity of transmitter output stages, and 
so on. The objective of this study is to obtain a computer-based prediction of noise and 
co-channel interference effects on the probability of bit error for M-ary Phase Shift 
Keying (MPSK) digital communication systems. To accomplish this, a model was 
developed using SIMULINK and the Communications Toolbox. The performance of the 
model has been verified by comparing the simulation results and the theoretical results for 
the bit error probability in the presence of additive white Gaussian noise (AWGN). The 
model has been subsequently modified to include co-channel interference, and the bit 
error probability of MPSK with AWGN and co-channel interference was obtained via 


simulation. 











Il. SIMULINK AND THE COMMUNICATIONS TOOLBOX 


SIMULINK is a program for modeling linear and non-linear dynamic systems in 


the time-domain. Models in SIMULINK are represented in block-diagram form and can 


- be assembled from block libraries and sub-libraries. Furthermore, SIMULINK allows the 


display of results and changes of certain model parameters without interrupting the 


simulation. Since SIMULINK is built upon the MATLAB numeric computation system, 


it offers direct access to the MATLAB workspace and MATLAB’s mathematical and 


engineering functions. SIMULINK’s main features are: 


modeling and analysis of dynamic systems, including linear, nonlinear, 
continuous, discrete, and hybrid 

flexible "open system" environment that allows addition of new blocks to 
Simulink 

seamless interface with MATLAB’s built-in math functions, 2-D and 3-D 
graphics, and add-on toolboxes for specialized applications 

choice of methods of running a simulation (menu-driven on-screen or batch- 
mode) 

an optimized computer platform implementation that ensures fast and accurate 
results 

unlimited model size. 


A typical SIMULINK session starts by either defining a new model or recalling a 


previously defined model and then proceeds to a simulation of the performance of that 
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mode]. In practice these two steps are often performed tiecstieke as the model designer 
creates and modifies a model to achieve the desired behavior [Ref.1]. 

The Communications Toolbox is a collection of MATLAB functions and 
SIMULINK blocks for model development and simulation in the communications area. 
The functions/blocks are organized in the following sub-categories: Data Source, Source 
Coding and Decoding, Error-Control Coding, Modulation and Demodulation, 


Transmission and Reception filters, Transmitting Channel, Multiple Access, and 


Synchronization and Utilities [Ref.1]. 








it. PHASE SHIFT KEYING 


Modulation of a sinusoidal carrier allows the conveyance of information in binary 
form. The baseband information bit stream modulates the carrier by means of discrete 
changes in the carrier amplitude (amplitude-shift keying or ASK) carrier frequency 
(frequency-shift keying or FSK) or carrier phase (phase-shift keying or PSK) or some 
combination of the above. Each different state of the carrier, known as a symbol, 


corresponds to one or more bits of baseband information [Ref 2.]. 


A. BINARY PHASE-SHIFT KEYING 


In binary phase-shift keying (BPSK), the phase of a constant amplitude carrier 
signal is switched between two values according to the two possible “messages” m and 
mz corresponding to binary 1 and Q, respectively. Normally, the two phases are separated 


by 180 degrees. If the sinusoidal carrier has amplitude A;, then the average energy per bit 


E, = = ArT, , and the transmitted BPSK signal can be represented as: 


S ppg (t) = a cos(277,t +0.) O<stsT, (binary 1) (1) | 
5 





S ppsg (t) = 2B cos(2a@f.t+m7+0,) O<t<T, (binary 0) (2) 
b 


= cos(2#ftt+O.) O<t<T, (binary0) © (3) 


b 





S ppsx (t) = — 











where E, is the average energy per bit, 7, is the bit duration, and a rectangular pulse 


shape p(t)=2[(¢ -7,,/2)/T,)] is assumed. For this signal set, the single waveform 


%,(t) = ie cos(27f.t) O<stsT, (4) 
b 


can be defined. Using this signal, the BPSK signal set can be represented as 


S apse = (VE, 0),-yE,4,0} | (5) 


B. SIGNAL CONSTELLATION FOR BPSK 





Figure 1. Signal Constellation for BPSK 
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A BPSK signal can be geometrically represented in “ signal space” as shown 
below. Such a representation is referred to as the symbol constellation and provides a 


graphical representation of the complex envelope of a BPSK signal for the two possible 











symbols. The x-axis of the constellation diagram represents the “in-phase” component of 
the complex envelope, and the y-axis represents the “quadrature” component of the 
complex envelope. The distance between the signals on a constellation diagram indicates 
how different the modulation waveforms are and determines how well a receiver can 
differentiate between all possible symbols when random noise is present. The larger the 
signal distance, the better the chance of correct symbol detection. 
C. PROBABILITY OF BIT ERROR FOR BPSK 

An important measure of performance used for digital modulation is the 
probability of symbol error, P, . It is often convenient to specify system performance by 
the probability of bit error P,,even when decisions are made on the basis of symbols 
rather than bits. The relationship between P, and P, for orthogonal signaling is: 


Pi, 


f= 6 
arr (6) 





For the BPSK modulation (M=2), the symbol error probability is equal to the bit 
error probability. When the signals are assumed equally likely and signal s,(t) G=1,2) is 
transmitted, the received signal, r(t), is equal to s,(t)+n(t), where n(t) is modeled as 
additive white Gaussian noise (AWGN). The antipodal signals (signals of equal 


amplitude and opposite polarity) s,(t) and s,(t) are (cf. Eq.5): 


s,(t) = JE, ¢,(t) O<tsT, (7) 
s,(t) = -JE, 9, (t) O<t<T, . (8) 





The decision stage of the detector will choose the s;(t) with the largest correlator output 
z,(t), or in this case of equal-energy antipodal signals, the detector, using the decision 
rule, decides: 

s@) if zM> 

s,(t) otherwise 
where ‘Yo denotes the decision threshold (equal to 0 for equally probable antipodal signals) 
and z (7) is the correlator output at time T. Two types of ietecuon error can - made. 
The first type of error takes place if s,(¢) is transmitted but the noise is such that the 
detector measures a negative value for z(T)and decides (incorrectly) that signal s,(¢) was 
sent. The second type of error takes place if signal s,(¢) is transmitted but the detector 
measures a positive value for 2(T) and decides (again incorrectly) that signal s, @) was 
sent. Therefore, the probability of bit error,P,, is the sum of two conditional 
probabilities: 

P, = P(s, \s,)P(s,)+ P(s, | s,)P(s,) | (9) 

_ For the case when the auibel probabilities are known , priori and the symbols are 
equally probable (which is mostly the case): 


P(s,) = P(s,) == | | (10) 


The expression for the bit error probability becomes: 


1 ; 
Fr ==P(s, I s,) rn: I's,) (11) 











Because of the symmetry of the probability density functions (pdf's) of the sum of signal 
and noise, P(z|s,) and P(z\.s,)the two conditional probabilities are also identical and 
P, = P(s, | s,) = P(s, | s,) | | (12) 
The probability of a BPSK bit error P, is numerically equal to the area under the “tail” of 
either pdf p(zls,) or p(zls,) that falls on the incorrect side of the threshold. We can 
therefore compute P, by either integrating p(z|s,) between the limits —°o and 7, or by 


integrating p(z|s,) between the limits 7, and co. Hence, 


P,= — { p(zls,)dz (13) 


19 =(Q +22) /2 
where the conditional pdf's p(zls,)@=12) are Gaussian with mean value a,, 


variance Oo, and the optimum threshold,7,, is (a4,+a,)/2. The probability of bit error 








for BPSK is: 
t 1 Z~Qy,3 
P; ae | —eewnl-5 ——+)? jdz 
¥,a+a,)/2 FoV 2H (14) 
If we introduce u=~——2 and 0,du=dz, 
| S, 

the integral simplifies to 

P, = f penile = Q(A—) (15) 

. w=(a-2,)!20) W 27 2 205 


where 0, is the standard deviation of the noise at the output of the correlator. The 


function Q(x) is defined as: 


{ exp(— du (16) 


Q(x) = I 


For equal energy antipodal signaling, such as the BPSK format in (5), the receiver output 
signal components are a, =a} E, when s,(t) is sent and a, = ~ JE, when s, (tf) is sent. 
For AWGN the noise variance o, at the correlator output is equal to N,/2 so that we 


can rewrite probability of bit error as follows [Ref.3]: 








r 1 u? 
P, = exp(-——)du 17 
, in exe (17) 
_ 2E, 
Pz= QO N, ) (18) 


D. QUADRATURE PHASE-SHIFT KEYING 


In quadrature phase-shift keying (QPSK), two bits are transmitted in a single 
modulation symbol. The phase of the carrier takes on one of four equally spaced values 


such as 0, 2/2, x, and 37/2, where each value of phase corresponds to a unique pair of 


message bits. A QPSK signal can be defined as 


Sopsx (t) = A cos[27f,t + G— va O<stsT. i=1,2,3,4. (19) 


where T,; denotes symbol duration (equal to twice the bit duration T,). By using 


trigonometric identities, Equation 19 can be written in the interval OS tS 7, as: | 


Sopsx (t) = = cos[(i —1) = cos(27f,t) — 1 = sin[ (i — v7 sin(27f.t) (20) 








LO 








If the basis functions ®,(t) = J2/ T, cos(27f,t) and ®, (t) = J2/ T, sin(27f.,t) are defined 
over the interval 0<¢<T. then the QPSK signal can be expressed in terms of the basis 


signals as [Ref.4]: 
Sorsx (t) = {of E, cosl(i- N10, (t) JE, sin[G- N50, ()} 1=1,23,4. (21) 


E. SIGNAL CONSTELLATION FOR QPSK 
The QPSK signal constellation is shown below. The constellation suggests that a 
QPSK signal can be thought of as a combination of two pairs of antipodal BPSK signals, 


with the two BPSK pairs orthogonal (in phase quadrature) to each other. 





Figure 2. Signal Constellation for QPSK 
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F. PROBABILITY OF BIT ERROR FOR QPSK | 
Since QPSK may be considered as a combination of two orthogonal BPSK 


signals, the bit QPSK error probability is identical to the BPSK bit error probability given 
by Equation 18. Note however that the QPSK symbol error probability is given [Ref.4] 


as: 


P= 20 pear Bs (22) 


G. SIMULINK MODEL OF MPSK DIGITAL COMMUNICATION SYSTEM 

In the simulation of communication systems, there are two options for 
_ representing the modulation/demodulation process: passband and baseband. In passband 
simulations, the carrier signal is included in the simulation model. By the Nyquist 
sampling theorem, the simulation sampling frequency must be at least twice the sum of 
the carrier signal frequency and the highest frequency component of the message in order 
to recover the message correctly. Since the frequency of the carrier is usually much 
greater than the highest frequency of the input message signal, a passband simulation 
requires a very large number of samples and excessive simulation times. The simulation 
time can be dramatically reduced by using baseband-equivalent simulation, also known 
as the low-pass equivalent method. A baseband equivalent model uses complex 
envelopes to reduce the simulation sampling rate and thus speed up the simulation. As 
the name suggests, the complex envelope is a complex function of time that carries 
information about both the time-varying amplitude and the time-varying phase of a 
passband signal. Since the amplitude of an MPSK signal does not vary with time, the 


complex envelope of an MPSK signal has a constant magnitude but a time-varying phase 
12 








angle of the complex envelope for the i-th symbol is 27i/M [Ref.1]. In this thesis the 


baseband equivalent model for MPSK, shown in Figure 3, was developed and applied to 


2PSK and 4PSK cases. 


In order to evaluate the effect of co-channel interference, the model in Figure 3 


was modified by adding another MPSK modulator and another message source. 


The MPSK model has the following main functional blocks: 
Sampled Read from Workspace | 
MPSK modulation and MPSK demodulation blocks 

AWGN channel 

K-step Delay 

Error Counter 


Triggered Write.to File. 


In addition to these blocks, auxiliary blocks such as multiplexer, oscilloscope, sum, pulse 


generator, switch, counter, and so on are also used. 


The Sampled Read From Workspace block reads a row of data for a workspace 


variable (a matrix) at every data sampling point. The output of this block is a vector 


whose length is equal to the number of columns of the workspace variable matrix. The 


possibility exists to offset (postpone) reading from workspace by a desired value such 
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MPSK mod 
baseband 








counter 


channel 


MPSK demod 
baseband 








Sampled read 
















from ae a ah 
Triggered wnite 
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Pulse 
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Figure 3. Model of MPSK Communication System 
that up to the user-defined offset time the output of this block is equal to a user-defined 


initial value. When a non-zero time offset is used, reading from workspace starts once 


the simulation time exceeds the offset. In the case that the last row of the workspace 


matrix is reached before completing the simulation (that is when the number of rows in 
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the workspace matrix is less than the number of simulation ies instants) the first row of 
the workspace variable may be read again if the user has selected the cyclic read option 
(otherwise the output is 0). 

The MPSK Mod block creates the complex envelope signal consisting of the real 
and imaginary parts of the signal. This block is a combination of the MPSK Map and CE 
(Complex Envelope) PM blocks. The digital input signals are in the range [0,M-1], where 
M is total number of different phases and the phase shift for input digit i is 2ni/M. The 
MPSK Mod output is a unit-magnitude complex analog signal whose ee values can 
be represented by M points on the unit circle in the complex plane. 

| The MPSK Demod block demodulates/demaps the complex input. This block is a 
combination of the MPSK Corr Demod and CE Min/Max Demap blocks. The MPSK 
Demod output is an integer in the range [0,M-1], where M is the number of phases. 

The AWGN (Additive White Gaussian Noise) Channel block adds white Gaussian 
noise with user-specified mean and alias to the signal being transmitted through this 
channel. This block uses the Gaussian White Noise Generator block as the noise source. 
The vector length (N) for the initial seed entry determines the output vector size. The 
mean and the variance of the noise are specified such that the mean can be either a vector 
of length equal to ie seed length (N) or a scalar, in which case all the elements of the 
noise vector share the same mean value. The covariance matrix can be one of the 
following: 

e AnN-by-N positive semi-definite matrix whose off-diagonal are the correlations 
e A length N vector, in which case the individual elements of the noise vector are 


uncorrelated but have unequal variances 
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e A scalar, in which case all the elements of the noise vector are uncorrelated but 
share the same variance. 

The K-Step Delay block delays its input (scalar) signal by K-sampling time 
intervals. This block function is equivalent to a z~ discrete time operator. 

The Error Counter block detects the difference between the port 1 signal and the 
port 2 signal at the times of the rising edge of the port 3 signal. If the difference between 
the port 1 and port 2 signals is larger than a user-defined tolerance, the error ‘counter 
increments its output by 1. The rising edge of the signal at port 4 resets the error counter 
back to 0. The Error Counter block accepts scalar signals only. 

The Triggered Write to Workspace block writes the input port 1 input variable to 
a workspace variable at the rising edge of the input port 2 trigger signal. There are six 
entries to this block: the workspace variable name, data type, the number of trigger pulses 
between saved data, the ceaanand row number, what to keep in the case of overflow, and 
the threshold value. The message signal to the input port 1 can be either a scalar or a 
vector while the trigger signal to the input port 2 must be a scalar. The saved workspace 
variable is a column vector if port 1 input is a scalar and a matrix if the port 1 input is a 
vector. In the later case the first element of the input signal vector is saved in the first 
column, the second element of the input signal vector is saved in the second column, and 
so on. The block saves the first record at the first rising edge of the trigger signal. The 
data can be saved either as strings or data. A pre-defined number limits the maximum 


row of the output variable (this number cannot be changed during simulation). _ 
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IV. BIT ERROR PROBABILITY CONVERGENCE 


In this chapter the convergence of the bit error probability for BPSK with 
AWGN is investigated. The probability of bit error obtained by simulation converges to 
the theoretical values as the “transmitted” number of symbols (and the simulation run 
time) increases. Convergence was tested for signal-to-noise (SNR) ratios with noise 
variances (the signal power was kept constant) of 0.5, 1, 2, 4, 9, and 16. The signal 
amplitude is denoted A (A=1), the symbol duration is denoted T (T= 1s), the sampling 
interval is denoted At (At=T/4=0.25s), and the noise variance is denoted o*. With this 


notation, the signal-to-noise ratio SNR in dB and the theoretical values of the bit error 


probability for BPSK are: 
AT 
SNR ap) = MOORS Ay? (23) 
2 SNR 
Ps=Q(V10 * ) (24) 


The theoretical values, calculated using Equations 23 and 24, are shown in Table I. 
Before showing the convergence results for each variance/SNR case, it is useful to explain 
the results for one ana For example, Table I shows the convergence results for the 
noise variance of unity. The table has the columns for: 

e the number of symbols used in the simulation 

e the theoretical value of the bit error probability 

e the bit error probability obtained by the simulation (for the corresponding number 


of symbols “transmitted”’) , and 
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e the difference between the theoretical and the simulation bit error probabilities. 
The table shows that the simulation result converges to the theoretical result as the 


number of “transmitted” symbols increases. 


Table I. Simulation Results for 67 =1 












The estimated bit error probability as a function of the symbol number for each of the 


simulations (from 10° to 10° symbols) is shown in Figures 4 through 7. 
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Figure 4. BER Convergence for 10° symbols 
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Figure 5. BER Convergence for 10* symbols 
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Figure 6. BER Convergence for 10° symbols 
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Figure 7. BER convergence for 10° symbols 
For the case of 1,000 symbols, the simulation gives the estimate of the probability 
of error as 0.025, as seen in Figure 4. Note that each time a simulation is repeated for the 
same number of symbols and the same variance but with a different seed value for the 
noise generator, a different result is in principle obtained. However, the average of all 
such simulation results becomes essentially a constant as the number of repeated 


simulations increases. When the number of transmitted symbols is increased to 10,000, a 
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different estimate for the bit error probability is obtained, 0.0237. If the number of | 
transmitted symbols is increased to 100,000, the bit error probability obtained is 0.0233 
and there is very little change past 3,000 or so symbols. Lastly, increasing the number of 
transmitted symbols to 1,000,000 yields the — accurate estimate for the bit error 
probability of 0.0231. 

Next, the number of transmitted symbols was selected as 10° and the simulation was run 
for the noise variances between 0.5 and 16. The results for the estimates of the bit error 


probability and their difference relative to the theoretical bit error probabilities are shown 


Theoretical P, Simulation P, Number of Difference 
Value Value Trials 


2.339*10° 


in Table II. 


0.0230 


0.1590 
0.2520 
0.3090 


Table I. Simulation Results for Various Variances 
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The convergence of the bit error estimates as functions of the numbers of transmitted 


symbols is shown in Figures 8 through 12. 
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Figure 8. BER Convergence for 0” 
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Figure 9. BER Convergence for 0” 
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Figure 12. BER Convergence for 0” = 16 
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The bit error probability estimates obtained by simulation for various SNR’s are shown in 


Figure 13 against the theoretical bit error probability for BPSK with AWGN. 
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Figure 13. BER vs. SNR 
As can be seen, the theoretical bit error probability and the Simulink results match each 


other closely. 
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V. SIMULATION VERIFICATION FOR BPSK AND QPSK WITH AWGN 


In this chapter the simulation results obtained from the (modified) Simulink model 
and a Matlab program will be presented and discussed. The Simulink model was modified 
to include co-channel interference (discussed in the next chapter). In this chapter, only the 
AWGN will be considered (the co-channel interference power was set to 0) for BPSK and 
QPSK modulations. The estimates for the bit error probability were obtained using both 
the Simulink model shown in Figure 14 and the Matlab program mpsk1.m hsted in the 
Appendix M. The mpsk1.m program uses Communications Toolbox functions rather than 
Simulink blocks to implement the simulations. Although a time-domain simulation of a 
communication system is in general possible using either Simulink block-diagrams or the 
equivalent Matlab functions, there is a very important difference between the two. The 
essential feature of Simulink is that simulations progress one step at the time, that is the 
input data is processed sequentially (one at a time), and are referred to as “point data” or 
“time-flow” simulations. On the other hand, Matlab simulations can process vectors 
(“blocks”) of data at each simulation step and are thus referred to as the “block-data” or 
“data-flow” simulations. Block-data simulations are substantially faster than point-data 
simulations. However, block-data simulations cannot be applied in the time-domain to — 
simulate feedback-systems since such systems require the responses to be calculated at 


each time instant before new inputs are applied. 
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Figure 14. Simulink MPSK Model with Co-channel Interference 
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For both BPSK and QPSK, only AWGN was considered and the co-channel interference 
was made negligible by setting the signal-to-interference ratio to +100 dB. Bit signal-to- 
noise (SNR) ratios in the range from —5 dB to +12 dB were used in the simulations. 
The simulation results were plotted using auxiliary Mathcad “scripts”: 

¢ plot_sim_noise2.mcd for Simulink BPSK simulations 

e plot_mat_noise2.med for Matlab BPSK simulations 

* plot_sim_noise4.mcd for Simulink QPSK simulations and 

e plot_mat_noise4.mcd for Matlab QPSK simulations. 
A. BPSK BIT ERROR PROBABILITY ESTIMATES 

1. BPSK Results for the Simulink Model 
The Simulink model was run from a Matlab script in order to automate the data 

input and data output operations for simulations run for multiple values of signal-to-noise 
ratio SNR. The input data file is created using Matlab program prep.m. Another Matlab 
program, sun_psk.m, reads the input data file, starts the simulation (executed in a loop for 
various values of SNR), and saves the simulation output data. The mput data used for the 
- simulation (and defined by running the prep.m Matlab script) are: 
Noise (1) or Noise and Interference (2): 1 
Symbol Duration: 1s 
Oversampling Factor: 2 
Minimum Bit Signal to Noise Ratio: -5 dB 
Maximum Bit Signal to Noise Ratio: 12 dB 


Number of values for SNR: 10 
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Minimum number of errors acceptable: 100 
Factor multiplying the error numbers: 2 
Maximum size of the random integer arrays: 10° 


File name to save data: SIMULINK_NOISE 2 


The output data saved in the ASCII format includes: 


a 


bit error probabilities (BER) 
signal-to-noise ratios (SNR) 
signal-to-interference ratios (SJR) (set to +100 dB for the noise-only case) 
numbers of errors observed 


numbers of transmitted symbols. 


The estimates (obtained from the Simulink model) for the BPSK bit error probability are 


shown in Figure 15 together with the theoretical BER versus the SNR . 
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Figure 15. BPSK BER versus SNR (Theory and Simulink Model Estimates) 

2. BPSK Results Obtained by Using Matlab Program MPSK1.M 

The same MPSK communication system that was modeled in Simulink using 
Simulink and Communications Toolbox blocks was implemented in Matlab using Matlab 
and Communications Toolbox functions. The same parameters were used for the 
simulation: 

Noise (1) or Noise and Interference (2): 1 

Symbol Duration: Is 

Oversampling Factor: 2 


Minimum Bit Signal to Noise Ratio: -5 dB 
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Maximum Bit Signal to Noise Ratio: 12 dB 

Number of values for SNR: 10 or 18 

Minimum number of errors acceptable: 100 

Factor multiplying the error numbers: 2 

Maximum size of the random integer arrays: 10° 

File name to save data: MATLAB_NOISE_2 
The simulation output data were saved and visualized using an auxiliary Mathcad script. 
The theoretical bit error probability and the estimates obtained from the Matlab ee 


are shown in Figure 16 as functions of the signal-to-noise ratio SNR. 
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Figure 16. BPSK BER versus SNR (Theory and Matlab Model Estimates) 
In this case, simulation results differ from theory. At closer inspection, one notices that 
the curve for the estimates is similar to the theoretical curve but shifted to the right. Such 
. a shift corresponds to an SNR “loss”. Indeed, if the SNR for the theoretical curve is 
reduced by 3 dB, the theoretical curve and the Matlab estimate match each other, as 
shown in Figure 17. Therefore, we may conclude that the Simulink model gives very 
accurate estimates for the bit error rate of BPSK with AWGN while the Matlab model 
experiences a 3 dB SNR loss (compensating for this SNR loss again produces accurate 
estimates). This seems to indicate a fundamental problem with the Communications 


Toolbox. 
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Figure 17. BPSK BER versus SNR (Theory for SNR-3dB and Matlab for SNR) 
B. QPSK BIT ERROR PROBABILITY ESTIMATES 
1. QPSK Results for the Simulink Model 
The Simulink model used for the BPSK is a general model that can be used for any 
MPSK. Therefore, the same model was used for QPSK by selecting the number of phases 
M (an input parameter) as four. The Simulink model results for the QPSK bit error 
probability estimates are shown in Figure 18, together with the theoretical bit error 


probability and its union bound given by [Ref.4] 


SNR,,; 
M log(M). ..——— . 4 
P. , .,(M,SNR,, »)=———Q(,| ve —"2.10 © sinc 25 
aoe sit _ dB) Mane: 10g(2) ve) (25) 
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as functions of the bit signal-to-noise ratio SNR»i. We note that the estimates differ 


somewhat from theory but that this difference diminishes as the bit SNR increases. 
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Figure 18. QPSK BER versus SNR (Theory and Simulink Model Estimates) 
2. QPSK Results Obtained by Using Matlab Program MPSK1.M 
The estimates for QPSK BER obtained using Matlab are shown in Figure 19, 
together with the ‘heretical bit error probability and its union-bound estimate given by 


Equation 25 as functions of the bit SNR. 
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Figure 19. QPSK BER versus SNR (Theory and Matlab Model Estimates) 
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VI. SIMULATION VERIFICATION FOR BPSK AND QPSK WITH 


ADDITIVE WHITE GAUSSIAN NOISE AND CO-CHANNEL INTERFERENCE 


A. BPSK BIT ERROR PROBABILITY ESTIMATES 

In the preceding chapter only additive white Gaussian noise (AWGN) was 
considered for BPSK and QPSK simulations. In this chapter both AWGN and co-channel 
interference are considered. A second MPSK Mod block fed by an independent message 
source generates the co-channel interference and is added to the sum of the signal and 
noise. In the Simulink model and the Matlab program a range of signal-to-noise ratios 
(SNR) and signal-to-interference ratios (SJR) are used to calculate the estimates of the bit 
error probability as functions of SNR and SJR for both BPSK and QPSK. The results are 
presented as families of curves, one with SNR as the variable and SJR as the variineti 
and the other with SJR as the variable as SNR as the parameter. 

The parameters for the BPSK simulations using both the Simulink model and 
Matlab were: _ 

Noise (1) or Noise and Interference (2): 2 

The Number of Phases: 2 

Symbol] Duration: Is 

Oversampling Factor: 2 

Minimum Bit Signal to Noise Ratio: -5 dB 

Maximum Bit Signal to Noise Ratio: 12 dB 


Minimum Bit Signal to Interference Ratio: -5 dB 
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Maximum Bit Signal to Interference Ratio: 12 dB 

Number of values for SNR: 10 

Number of values for SJR: 10 

Minimum number of errors acceptable: 100 

Factor multiplying the number of error: 2 

Maximum size of the random integer arrays: 10° 

File name to save data: NOISE&INTERFERENCE_2 

The probabilities of bit error were saved as a BER matrix whose rows correspond 
to constant values of SJR and whose columns correspond to constant values of SNR. The 
selected ranges for SNR and SJR were from —5 dB to +12 dB. For each bit error 
probability estimate (each matrix entry) at least 100 bit errors were observed. The 
simulations ran until the error counter exceeded the minimum specified ambero error, in 
this case 100. Upon exceeding the minimum number of required bit errors, the simulation 
was restarted for another combination of SNR and SJR (to create another BER matrix 
entry). High values of SNR or SJR correspond to low noise and interference powers, 
which in turn implies low probabilities of bit error. Therefore, the time required for the 
simulation to produce at least 100 errors increases as the row and column indices (SJR 
and SNR) increase, and the BER matrix entries in general decrease with increasing the 
row and column indices. 

For QPSK simulations the same parameter values were used as for BPSK except 
for the number of phases M, which are _ for QPSK. The BPSK simulation results were 


plotted using auxiliary Mathcad programs plot_sim_noise&int2.mcd for the Simulink 
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model and plot_mat_noise&int2.mcd for Matlab program results. Similarly, the QPSK 
simulation results were plotted using auxiliary Mathcad programs 
plot_sim_noise&int4.mcd for Simulink model and plot_mat_noise&int4.mcd for the 
Matlab program. 

1. BPSK Results for the Simulink Model 

The results for the bit error probability of BPSK with noise and siaianel 
interference obtained using the Simulink model are shown in Figure 20 (BER versus SNR 
with SJR as parameter) and in Figure 21 (BER versus SJR with SNR as parameter). The 


increments for both SNR and SJR are 1.889 dB, starting from —5 dB. 
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Figure 20. BPSK BER versus SNR with SJR as Parameter (Simulink Model Estimates) 

The “lowest” curve represents the theoretical bit error probability for BPSK with AWGN 
only. The increase of BER with decreasing SJR is evident from Figure 20. The BPSK 
probability of bit error as a function of SJR, with SNR as a parameter starting at —5 dB 


and with 1.889-dB increments, is shown in Figure 21. 
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Figure 21. BPSK BER versus SJR with SNR as Parameter (Simulink Model Estimates) 
The “step function” represents the noise-free case, that is the case of co-channel 
interference only. In such a case the probability of error is either 0.5 when the co-channel 


interference power is larger than the signal power (negative SJR) or 0 when the signal 
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power is larger than the co-channel interference power. (In order to use the logarithmic 
scale, the value of 0 for positive SJR has been replaced by a small non-zero value). 

2. BPSK Results Obtained by Using Matlab Program MPSK1.M 

The results for the bit error probability of BPSK with noise and co-channel 
interference obtained using Matlab program are shown in Figures 22 and 23 (BER versus 
SNR with SJR as parameter starting at —-5 dB and with 1.889 dB increments) and m 
Figure 24 (BER versus SJR with SNR as parameter starting at -5 dB and with 1.889 dB 
increments). Almost 3 dB difference can be noticed between simulation results and the 


theoretical value. This difference was explained in the preceding section. 
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Figure 22. BPSK BER versus SNR with SJR as Parameter (Matlab Model Estimates) 


Figure 23 shows the BPSK probability of bit error as a function of SJR, with SNR as a 
parameter. Again, the difference between SNR curves is 1.889 dB as in SJR case. The 
“step function” again represents the noise-free case, that is the case of co-channel 


interference only. 
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Figure 23. BPSK BER versus SJR with SNR as Parameter (Matlab Model Estimates) 
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B. QPSK BIT ERROR PROBABILITY ESTIMATES 

1. QPSK Results for the Simulink Model 

The — for the bit error probability of QPSK with noise and co-channel 
interference obtained using the Simulink model are shown in Figure 24 (BER versus SNR 
with SJR as parameter starting at —5 dB and with 1.889 dB increments) and in Figure 25 
(BER versus SJR with SNR as parameter starting at —5 dB and with 1.889-dB 


increments). 
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Figure 24. QPSK BER versus SNR with SJR as Parameter (Simulink Model Estimates) 
The “lowest” curve represents the theoretical bit error probability for QPSK with AWGN 
only. The increase of BER with decreasing SJR (for positive SNR) is evident from Figure 


24. 
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The QPSK probability of bit error as a function of SJR, with SNR as a parameter 
starting at —5 dB and with 1.889-dB increments, is shown in Figure 25. The “step 
function” again represents the noise-free case, that is the case of co-channel interference 


only. 
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Figure 25. QPSK BER versus SJR with SNR as Parameter (Simulink Model Estimates) 
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As can be seen in Figure 25, the curves for the bit error probability tend to the noise-free 
step function as the SNR increases. Also, the curves for the negative SJR tend to 0.5, 
meaning that the bit error probability is dominated by the interference for negative SJR, as 
expected. 

2. QPSK Results Obtained by Using Matlab Program MPSK1.M 

The estimates for QPSK BER obtained using Matlab are shown in Figure 26 
together with the theoretical bit error probability as a function of the bit SNR and with 
SJR as parameter starting at -5 dB and with 1.889 dB increments. Comparing the BER 


curves in Figures 24 and 26, we note that they differ slightly. 
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Figure 26. QPSK BER versus SNR with SJR as Parameter (Matlab Model Estimates) 

The QPSK probability of bit error as a function of SJR, with SNR as a parameter 
starting at —5 dB and with 1.889-dB increments, is shown in Figure 27. The “step 
function” again represents the noise-free case, and the curves are seen to tend to the step 


function as the SNR increases. 
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Figure 27. QPSK BER versus SJR with SNR as Parameter (Matlab Model Estimates) 


From these figures, we see that co-channel interference acts like noise. When co-channel 


interference is high (SJR is low), the bit error probability is high and when it is low (SJR is 


high) the bit error probability is low provided SNR is also high. 
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From an examination of Figures 26 and 27, we see that we can maintain je 
for SIR< 12 dB when SNR28.22 dB. Up to a point, further increasing SNR allows SJR 
to further decrease. A similar result is observed for BPSK, as can be seen from an 


examination of Figures 20 through 23. Hence, the importance of AWGN on the effect of 


co-channel interference cannot be overlooked. 


=p) 
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VIL CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The objective of this thesis was to model and simulate MPSK communication 
systems in the presence of additive white Gaussian noise and co-channel interference. The 
models were implemented as Simulink block-diagram models, using Simulink and 
Communications Toolbox blocks, as well as Matlab programs, using Matlab and 
specialized Communications Toolbox functions. Baseband equivalents of passband MPSK 
systems were modeled using the complex envelope concept. The results for the bit error 
probability were obtamed for BPSK (2PSK) and OPSK (4PSK) modulation types with 
additive white Gaussian noise only and with the sum of additive white Gaussian noise and 
co-channel interference. Furthermore, the convergence (with _— numbers of 
“transmitted” symbols) of bit error probability estimates obtained by simulation to the 
theoretical bit error probabilities has been verified for several cases of additive white 
Gaussian noise with different variances. 

The simulation results For AWGN only are in eee very close to the theoretical | 
values and the estimates of the bit error rates converge to the theoretical values as the 
numbers of transmitted symbols increase. However, it has been observed that the Matlab 
model requires a 3 dB SNR correction for BPSK in noise case only. With this SNR 
correction, the simulation results for the bit error probability match the theoretical results 


and the results for the two implementations (SUMULINK and MATLAB) match each 


i 











other. The BPSK and QPSK are both sensitive to co-channel interference almost to the 
same extent. 

As a general rule, the sensitivity of both BPSK and QPSK to co-channel interference 
can be reduced by increasing the SNR. Hence, the importance of AWGN on the effect of 
co-channel interference cannot be overlooked. Of course, if the SJR becomes too small, 
then performance is poor regardless of the SNR. 

B. RECOMMENDATIONS 

This research can be continued to include forward error correction and adjacent 
channel interference in the developed models. Since only linear MPSK systems were 
| considered, a future study may address transmitter non-linearity. Theoretical derivation of 
bit error probability can be obtamed for MPSK with AWGN and interference which can be 
either phase locked to the signal or with some given distribution (such as Uniform or 


Gaussian) of the phase difference between the interference and the signal. 
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APPENDIX A. BATCH_RUN.M MATLAB PROGRAM FOR CONVERGENCE 


TEST OF SIMULATION USING BPSK 


% This program runs multiple simulations for convergence of simulink’ 


clear 


num_symbols = input(’Enter the number of symbols [1000]: ‘’); 
if isempty(num_symbols), num_symbols = 1000; end 


tic 
sim(’mpsk7_z’,num_symbols) ; 
disp(’VAR 0.5 Done!’) 

toc 


tic 
sim(’mosk7_1’,num_symbols) ; 
disp(’VAR 1 Done!’) 

toc 


tic 
sim(‘’mosk7_2’,num_symbols) ; 
disp(’VAR 2 Done!’) 

toc. 


cL | 
sim(‘mpsk7_4’,num_symbols) ; 
disp(’VAR 4 Done!’) 

toc 


Lic 
sim(’mpsk7_9’,num_symbols) ; 
disp(’VAR 9 Done! ’) 

toc 


alk e: 
sim(’mosk7_16’,num_symbols) ; 


disp(’VAR 16 Done!’) 
toc 
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APPENDIX B. BER_PLOT.M MATLAB PROGRAM FOR PLOTTING THE 


FIGURES OF THE SIMULATION FOR EACH VARIANCES 


% This loads two files and plots the result 
clear 

num bits = input(’Enter the number of bits [1000]: '); 
if isempty(num_bits), num_bits = 1000; end 
load d:\errnumi.dat 

load d:\bitnuml.dat 

plot(bitnuml, errnuml ./ bitnuml, ‘r-’) 

grid 

xlabel (’Bit Number’ ) 

ylabel(’Bit Error Probability’) 

title(’Bit Error Probability for variance 1’) 
[num_errors dummy] = size(errnuml) ; 


disp(’noise variance=1’) 
BER = num_errors / num_bits 


6] 
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APPENDIX C. PLOT_SIM_NOISE2Z MATHCAD PROGRAM FOR RESULTS OF 


NOISE ONLY CASE BPSK USING SIMULINK MODEL 


This program gives the probability of bit error vs. signal to noise ratio using outputs obtained 
from SIMULATION test for noise only BPSK. 


BER := READPRN ( "“h:\veys\simnoise2.ber” _—) 


SNR := READPRN ("“h:\veys\simnoise2.snr" —~) 


_ T 
SNR = SNR m :=0.. rows (SNR) — 1 
| SNR gp 
| _ 1 I 10 ies 
P @ bit theory (SNR dB) aa 1 — erf| —-4/2-10 | SNR gp = -5,-4.9.. 13 
,0.238095 , 





0.01 


1°10 > 


1*10 + 


"10 > 





1910 © 


1%10 7 


i%10 ° 


1910°” 





-10 om 
u! 332931 "10 2 1*10 10 
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APPENDIX D. PLOT_MAT_NOISE2 MATHCAD PROGRAM FOR RESULTS 


OF NOISE ONLY CASE BPSK USING MATLAB PROGRAM 


This program gives the probability of bit error vs. signal to noise ratio using outputs obtained from 
MATLAB test for noise only BPSK. 


BER '= READPRN("h:\veys\matnoise2.ber" ) 


SNR ‘= READPRN("h:\veys\matnoise2.snr"_——) 


SNR != SNR! 


m:=0.. rows(SNR) — 1 





BER 11 0° 








o-oo 
P « bit theory(SNR ap) wat 











plot with 3 db error 
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plot 3 dB corrected 
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APPENDIX E. PLOT_SIM_NOISE4 MATHCAD PROGRAM FOR RESULTS OF 


NOISE ONLY CASE QPSK USING SIMULINK MODEL 


This program gives the probability of bit error vs. signal to noise ratio using outputs obtained 
from SIMULATION test for noise only QPSK. 


Mis4 
BER :=READPRN ("h:\veys\simnoise4.ber" —_) 
SNR '= READPRN ("h:\veys\simnoise4.snr" —) 


SNR := snr! m :=0.. rows (SNR) — 1 


















BER :=t™M_ per SNR bit_aB_sim | “SNR - > 
2M-— 
Q(x) =4-[1- erf =] 
q2 
SNR pit dB 
Q¢(M.SNRpit ap) =Q 6 wold) 1g : 


M-—1 log(2) 





Pe bit QAM (M.SNR bit_dB) *= 
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E o¢(O.M,SNR pit) *=1+ erf(p(0,M, SNRpjt)) 


2 
© (0.M, SNR) =9(8.M,SNR yy) <P (8M, SNR 
SNR pit 
~10 10 log(M) 
log( 2) 
p 9(8.M,SNR,j,) =2—_—__—- 


2a 


Zee (0.m. sum) 


Ax 








SNR pit dB (=-9,-6.9.. 10 
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0.410547, 
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APPENDIX F. PLOT_MAT_NOISE4 MATHCAD PROGRAM FOR RESULTS 


OF NOISE ONLY CASE QPSK USING MATLAB PROGRAM 


This program gives the probability of bit error vs. signal to noise ratio using outputs obtained from 
MATLAB test for noise only QPSK. 


ea Me a Se 
oe ty 

mY, oe 

+4 i 

Re a ae EE 


BER:=READPRN "h:\veys\matnoise4.ber" ) 


SNR ‘=READPRN "h:\veys\matnoise4.snr"_ ) 
SNR:=SNR! 


m:=Q.. rows(SNR) — 4 


vet 








6 log) 6 


EE: 


M-—1  log(2) 





Pe bit. OAM (M.SNRpit ap) = — Q(M,SNRpit ap) I" 1-7 or sn) 
ce AM 


AM 
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APPENDIX G. PLOT_SIM_NOISE&INT2Z MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE BPSK USING SIMULINK 


MODEL 


This program gives the probability of bit error vs. signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from SIMULATION test for BPSK. 


BER !=READPRN ("“h:\veys\simnoint2.ber" +) 
SNR :=READPRN ("“h:\veys\simnoint2.snr"_—s) 


SJR :=READPRN ("“h:\veys\simnoint2.sjr"—) 


SNR !=SNR! sR :=SyR! 
m :=0.. rows (SNR) — 1 
n :=0.. rows (SJR) — 1 
| SNR gp 
od I 10 
P theo (SNR gp) =] = ef 242-10 


SNR yp *=-5.. 12 
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e-© SJR=-0.67 dB 
mene SJR=2.56 dB 
—— SJR=4.44 dB 
wweme SIR=6.33 dB 
—— SJR=8.22 dB 


eps =10> 
BERint only (SIR dB> eps ) = if( SJR dp<0-0.5, eps ) 


SIR gp 1=-5,- 4.9.12 
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APPENDIX H. PLOT_MAT_NOISE&INT2 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE BPSK USING MATLAB 
PROGRAM 


This program gives the probability of bit error vs. Signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from MATLAB test for BPSK. 


BER :‘=READPRN ("h:\veys\matnoint2.ber" _—+) 
SNR :=READPRN ("h:\veys\matnoint2.snr”_—+*) 


SJR ‘=READPRN ("h:\veys\matnoint2.sjr"—_—+) 





SNR !=SNR- SIR :=SyRT 
m :=0.. rows (SNR) — 1 
n :=Q.. rows (SIR) — 1 
SNR gp 
P theo (SNR gp) =—]1-erf] L210 © ‘= 
theo (SNR ap ) =>] 1— erf| 2 SNR gp !=-5..12 
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eps :=10°> 


BER (SIR gp, eps) #=if( SIR gp <0,0.5, eps ) 
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SJR gp =-5,-4.9.. 12 



























0,n . 
XK es ONS, ee a | 
BER, nm a 
BHo Se 
BER, ene Cenies Eee a pee 
if nS SO a = 
BER; 
o-o 2 
BER) in 0.0) eee OO 
‘ eae NE EEE eat nS A, CL I, RR ENT 
Gace : 
ee 
BERS on fare ene Ne eee ee 
BE aaa (EAI) SR, ET Sa, 
‘BER. ho ~~ 
eee. —3 | | 
ee 1%10 Rete Se pe pees 5 pe ee Ne | 
fe A ee eR ee Ete oe) 
BERS _ Noms 
BENS a 





—— SNR=INFINITY 
%* SNR=-5.00 dB 
ES SNR=-3.11 dB 
—2~- SNR=-1.22 dB 
@© SNR=0.67 dB 
eee SNR=2.56 dB 
wnemeee §=SONR=-4.44 dB 
~~ SNR=6.33 dB 
— SnR=8224B 
------- SNR=12.00 dB 


81 





82 











APPENDIX I. PLOT_SIM_NOISE&INT4 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE QPSK USING SIMULINK 


MODEL 


This program gives the probability of bit error vs. signal to noise ratio and probability of bit 
error vs.signal to jamming ratio using outputs obtained from SIMULATION test for QPSK. 


M := 

BER :=READPRN ("h:\veys\simnoint4.ber" —) 
SNR :=READPRN ("“h:\veys\simnoint4.snr" Ss) 
SJR ‘=READPRN ("“h:\veys\simnoint4.sjr" ss) 





SNR :=SNR! —3 SIR :=SIR BER == — -BER 
m :=0.. rows (SNR) — 1 


n :=0.. rows (SJR) — 1 


Q(x) e2-f1— ert{ 
2 
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eps :=10°> 


BERint only (SIR gp. eps ) =i (SIR dp <9, 9.5, eps | : 


SIR gp =-5,-4.9.. 12 
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APPENDIX J. PLOT_MAT_NOISE&INT4 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE QPSK USING MATLAB 
PROGRAM 


This program gives the probability of bit error vs. signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from MATLAB test for QPSK. 


M :=4 
BER :‘=READPRN ("h:\veys\matnoint4.ber" ——) 
SNR :=READPRN ("h:\veys\matnoint4.snr" —s) 


SJR :=READPRN ("h:\veys\matnoint4.sjr" +) 


SNR :=SNR! syR:=SJR! ‘BER !=BER 
m :=0.. rows (SNR) — 1 Q(x) mi 1— ert|_* 
n :=0.. rows (SJR) — 1 2 2 
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APPENDIX K. PREPAPE.M MATLAB PROGRAM TO ENTER THE SYSTEM 
PARAMETERS TO BE USED BY THE SIMUINK MODEL 


%%% This prepares the data file for simulink runs 
clear 


noise_only = menu(’ Select: ’... 
” NOISE’... 
’ NOISE and INTERFERENCE ; 


num_levels = input(Enter the number of frequencies M [2]: ); 
if isempty(num_levels), num_levels = 2; end 


T_sym | =input(Enter the symbol duration T [1]: ’); 
if isempty(T_sym), T_sym = 1; end 


oversampling = input(Enter the oversampling factor [2]: ’); 
if isempty(oversampling), oversampling = 2; end 


min_SNR = input(Enter the MIN Signal to Noise ratio[-5 dB]: ); 
if isempty(min_SNR), min_SNR = -5; end 


max_SNR = input(Enter the MAX Signal to Noise ratio [12 dB]: ); 
if isempty(max_SNR), max_SNR = 12; end 
if noise_only ~= 1 
min_SJR = input(Enter the MIN Signal to Interference ratio [-5 dB]: ); 
if isemptyGnin_SJR), min_SJR = -5; end 
max_SJR = input(Enter the MAX Signal to Interference ratio[12 dB]: )); 
if isempty(max_SJR), max_SJR = 12; end 
end 
if min_SNR == max_SNR 
num_noise = 1; 
else 
num_noise = input(Enter the number of values for SNR [10]: ); 
if isempty(num_noise), num_noise = 10; end 
end 


if noise_only ~= 1 
if min SJR == max_SJR 
num_jam = 1; | 
else 
num_jam = input(Enter the number of values for SJR [10]: 9; 
if isempty(num_jam), num_jam = 10; end 
end 
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else 
num_jam = 1; 
end 


min_errors = input(Enter the min number of errors acceptable [100]: ); 
if isempty(min_errors), min_errors = 100; end 


error_factor = input(Enter the factor multiplying the number of errors [2]: 9; 
if isempty(error_factor), error_factor = 2; end | 


initial_num_symbols = error_factor*min_errors; 


max_randint = input(Enter the maximum size of the random integer arrays [10“6]: ); 
if isempty(max_randint), max_randint = 10*6; end 


file_name = input(Enter the file name to save data [no ext]: 4's); _ 


save sun_data 
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APPENDIX L. MPSK.M MATLAB PROGRAM TO START SIMULATION | 


So%% This runs MPSK co-channel interference with additive noise 
clear 


load sun_data %This loads all the input data 


delta_t =T_sym/(2*oversampling) 

seeds =randint(3,1,1000); 

signal_seed = seeds(1); 

noise_seed = seeds(2); 

interf_seed = seeds(3); 

initial_num_symbols = error_factor*min_errors; 


tic 


if num_noise > 1 

delta_SNR = (max_SNR - min_SNR) / (num_noise -1); 
else 

delta_SNR = 0; 
end | 
SNR = min_SNR + [0:num_noise - 1]*delta_SNR; | 
Jonoise_var_vect = T_sym/(2*delta_t) .* 10 .* (-SNR/10); - 
noise_var_vect = 10 .4 (-SNR/10); 


if noise_only ~= 1 
if num_jam > 1 
delta_SJR = (max_SJR - min_SJR) / (num_jam -1); 
else 
delta_SJR = 0; 
end 
SJR = min_SJR + [O:num_jam - 1]*delta_SJR; 
interf_gain_vect = 10 .* (-SJR/20); 


else 
SJR = - 100; % There is no Jamming so the SJR in dB is -infinity 
interf_gain_vect = 0; 

end 


BER = zeros(num_noise,num_jam); 
total_symbols = initial_num_symbols; — 
for noise_case = 1:num_noise 
noise_var = noise_var_vect(noise_case); 
for jam_case = l:num_jam 
enough_errors = 0; 
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num_ err = 0; 

num_symbols = initial_num_symbols; 
rand_int =min({num_symbols max_randint]); 
total_symbols = num_symbols 


while enough_errors ~= | 
interf_gain = interf_gain_vect(jam_case); 
rand_int=min({10*6 num_symbols]); 
clear error_number 
sim(’pskco_bm’,num_symbols) 
[new_errors err_cols] = size(error_number); 
num_err = num_err + new_errors; 


if num_err == 
num_symbols =num_symbols*min_errors; 
and_int =min([num_symbols max_randint]); 


total_symbols = total_symbols + num_symbols; 

elseif (num_err > 0 & num_err < min_errors) 

num_symbols = (min_errors - num_err)*ceil(total_symbols/num_err); 
total_symbols = total_symbols + num_symbols; 

else 
enough_errors = 1; 
end 
end 


number_of_errors(noise_case,jam_case) = num_err 
number_of_symbols(noise_case,jam_case) = total_symbols 
BER(noise_case,jam_case) = num_err / total_symbols 

end 

eval([’ save ’ file_name ’.ssnr SNR -ascii]); 

eval([’ save ’ file_name *.sjr SJR -asci1]); 

eval([’ save ’ file_name ’.ber BER -ascii’)}); 

eval([’ save ’ file_name ’ner number_of_errors -ascii’]); 

eval([’ save ’ file_name ’.nsy number_of_symbols -ascii]) 


end 
toc 
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APPENDIX M. MPSK1.M MATLAB PROGRAM FOR USING MATLAB 


PROGRAM OPTION VICE SIMULINK MODEL 


JoJo Vo Te Wo Ve To So Yo Vo Ve Yo Vo Po Vo Ue Vo Ve Yo To Ne To Vo Yo Vo Po Vo Vo Ve 


%%% Matlab Version for MPSK GV 
%I% Dr. Jovan Lebaric abe lee o 
%%% June 14, 1998 Vo Vo 


Wo Vo Io Vo Vo Vo Vo Vo To To Vo Vo Go Fo Vo Go Vo Lo Vo Po Vo Vo LoL PRU 
clear 


%%KI% User Input %%% 
To Vo Toe Yo Ue Le Vo Vo Vo Vo Lo Vo Yo Vo Vo Vo Vo 


j= sqrt(-1); 


noise_only = menu(’ Select: ’... 

” NOISE’... 

>NOISE and INTERFERENCE };M = input(Enter the number of phases 
(power of 2) M [2]: 9; 
if isempty(M), M = 2; end 


T_sym = input(Enter the symbol duration T [1]: ); 
if isempty(T_sym), T_sym = 1; end 


oversampling = input(Enter the oversampling factor [2]: }); 
if isempty(oversampling), oversampling = 2; end 


min_SNR_bit = input(Enter the Minimum BIT Signal to Noise ratio in dB[-5 dB]: ); 
if isempty(min_SNR_bit), min_SNR_bit = -5; end 


max_SNR_bit = input(Enter the Maximum BIT Signal to Noise ratio in dB [+12 dB]: ); 
if isempty(max_SNR_bit), max_SNR_bit = 12; end 


if noise_only ~= 1 
min_SJR_bit = input(Enter the Minimum BIT Signal to Interference ratio in dB [-5 
dB]: ’; 
if isempty(min_SJR_bit), min_SJR_bit = -5; end 
max_SJR_bit = input(Enter the Maximum BIT Signal to Interference ratio in dB 
[+12 dB]: ); | 
if isempty(max_SJR_bit), max_SJR_bit = 12; end 
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end 


if max_SNR_bit == min_SNR_ bit 
num_noise = 1; 
else 
num_noise = input(Enter the number of values for SNR [18]: ); 


if isempty(num_noise), num_noise = 10; end 
end 


if noise_only ~= I 
if min_SJR_bit == max_SJR_bit 
num_jam = 1; 
else 
num_jam = input(Enter the number of values for SJR [10]: }); 
if isempty(num_jam), num_jam = 10; end 
end 
else 
num_jam = 1; 
end 


min_errors = input(Enter the min number of errors acceptable [100]: i 
if isempty(min_errors), min_errors = 100; end 


error_factor = input(Enter the factor multiplying the number of errors [2]: ); 
if isempty(error_factor), error_factor = 2; end 


initial_num_symbols = error_factor*min_errors; 


max_randint = input(Enter the maximum size of the random integer arrays [106]: ); 
if isempty(max_randint), max_randint = 106; end 


file_name = input(Enter the file name to save data [no ext]: ’,’s); 


%Ie% Parameters for MPSK %2%% 


%Ve% Symbol Frequency %%% 
Fd = 1/T_sym; 


%%% Sampling Frequency %%% 
Fs = oversampling* Fd; 


%%% Sampling Interval 2%% 
delta_t = 1/Fs; 
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%%% Seeds for Random Integers %%% 
seeds = randint(2,1,1000); 
signal_seed = seeds(1); 

interf_seed = seeds(2); 


initial_num_symbols = error_factor*min_errors; 
tic 


%IWeG% SNR Array 6% 
if num_noise > 1 

delta_SNR_bit = (max_SNR_bit - min_SNR_bit) / num_noise -1), 
else 

delta_SNR_bit = 0; 
end 
SNR_bit = min_SNR_bit + [O:num_noise - 1]*delta_SNR_bit;’ 
SNR_symbol = SNR_bit + 10*log10(log(M)/log(2)); 


%%% Noise Variance Array 26% 
noise_var_vect = 10.“ (-SNR_symbol/10); 
%noise_var_vect = T_sym/ (2*delta_t) * 10 .* (-SNR_symbol/10); 


%%e% SIR Array CHR 
if noise_only ~= 1 


if num_jam > 1 

delta_SJR_bit = (max_SJR_bit - min_SJR_bit) / (num_jam -1); 
else 

delta_SJR_bit = 0; 
end 


SJR_bit = min_SJR_bit + [O:num_jam - 1]*delta_SJR_bit; 
SJR_symbol = SJR_bit + 10*log10(log(M)/log(2)); 
interf_gain_vect = 10 .* (-SJR_symbol/20); 


else 


SJR_symbol = 100; | 
SJR_bit =100; % There is no Jamming so the SJR is large (infinity) 


end 


BER = zeros(num_noise,num_jam); 
total_symbols = initial_num_symbols; 
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Teo To Te Fo Vo Ve To Fo Fo Wo Yo Yo Vo To Fo Vo Yo Yo We Fo Vo Uo Vo To Vo Ge To Fo Ve Lo Le Vo Vo Vo Vo Ge 
%I%e% Noise and Interference Loops %%% 


for noise_case = 1:num_noise 
noise_var = noise_var_vect(noise_case); 
for jam_case = l:num_jam 


enough_errors = 0; 


num_err =Q; 
num_symbols = initial_num_symbols; 
rand_int = min([num_symbols max_randint}); 


total_symbols = num_symbols 


while enough_errors ~= 1 


rand_int = min([10“6 num_symbols]); 
input_sequence = randint(rand_int, 1,[0 M-1], signal_seed); 


%%% Modulation %%%% 
mod_MPSK = dmodce(input_sequence, Fd, Fs, ‘psk’, M); 


% Io %o Mo Vo Yo To Yo To Yo Yo Vo Vo Wo Yo Yo Ue Vo To Yo Vo We Yo To To To Vo Yo Vo Ve Veo 

%%% Additive Gaussian Noise %%% 

randn(’state’ssum(100*clock)); | 

complex_noise = randn(rand_int*oversampling, 1)*sqrt(noise_var); 

Serandn(‘state’,ssum(200*clock)); 

%complex_noise = complex_noise + 
j*randn(rand_int*oversampling, 1)*sqrt(noise_var); 


mod_MPSK = mod_MPSK + complex_noise; 


%%KI% Additive Interference %%% 
if noise_only ~= 1 
interf_gain = interf_gain_vect(jam_case); 
mod_MPSK = mod_MPSK + interf_gain*randint(rand_int*oversampling, 1,[0 
M-1], interf_seed); 
end - 


Teo To We Yo Ve Vo Yo Vo Vo Ve Lo Ve Lo Go Vo Ke Ve Vo Vo 
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%K% Demodulation %%% 
demod_MPSK = ddemodce(mod_MPSK, Fd,Fs, ‘psk’,M); 


ToT Vo% Vo % Le Io Go Ve Ve Vo Go % 
%I%% Errors GG 


[symbols_in_error dummy_col] = find(demod_MPSK - input_sequence),; 
[new_errors dummy_col] = size(symbols_in_error); 
num_err = num_err + new_errors; 


if num_err == 
num_symbols = num_symbols*min_errors; 
rand_int =min({num_symbols max_randint]); | 


total_symbols = total_symbols + num_symbols; 
elseif (num_err > 0 & num_err < min_errors) 
num_symbols = (min_errors - num_err)*ceil(total_symbols/num_err); 
total_symbols = total_symbols + num_symbols; 
else | 
enough_errors = 1; 
end 


end 
number_of_errors(noise_case,jam_case) = num_err 
number_of_symbols(noise_case,jam_case) = total_symbols 
BER(noise_case,jam_case) = 0.5 * M/(M-1) * num_err/ total_symbols 

end 

eval([’ save ’ file_name ’snr SNR_bit -ascii7]); 

eval([’ save ’ file_name ’.sjr SJR_bit -ascii’]); 

eval([’ save ’ file_name ’.ber BER -ascii’)]); 

eval({’ save ’file_name ’.ner number_of_errors -ascii’]); 

eval([’ save ’ file_name ’nsy number_of_symbols -ascii); 


end 


toc 


99 


100 





LIST OF REFERENCES 


. Communications Toolbox for Use with MATLAB and SIMULINK, The 








MATHWORKS Inc. 


. Hwei P.Hsu, Schaum’s Outlines, Analog and Digital Communications, McGraw-Hill, 
1993. | 


. Bernard Sklar, Digital Communications Fundamental and Applications, Prentice Hall, 
New Jersey, 1988. 


_ John Proakis, Digital Communications, McGraw-Hill, New York, 1995. 


101 





102 








BIBLIOGRAPHY 


. Bernard Sklar, Digital Communications Fundamental and Applications, Prentice Hall, 
New Jersey, 1988. 


. John Proakis, Digital Communications, McGraw-Hill, New York, 1995. 
. LA. Glover & P. M.Grant, Digital Communications, Prentice Hall Europe, 1998. 


. Communications Toolbox for Use with MATLAB and SIMULINK, The 
MATHWORKS Inc. 


. HweiP.Hsu, Schaum’s Outlines, Analog and Digital Communications, McGraw-Hill, 
1993. 


. Dr. Jovan Lebaric, Notes for EC4550 (Digital Communications), Naval Postgraduate 
School, 1998 (unpublished). 


103 


104 











10. 





INITIAL DISTRIBUTION LIST 


Defense Technical Information Center. ...................cccececececececenceccces 


8725 John J. Kingman Rd., STE 0944 
Fort Belvoir, VA 22060-6218 


Dyciey OX Tar y see sses ect atpcieseses ceases ee pee 


Naval Postgraduate School 
411 Dyer Rd. 
Monterey, CA 93943-5101 


GAP Vig tetsu Mal O66 (ol 2) Ce aa naat a tea ee sea oon Dy BPE etre OE oENn sree OORT See pete 


Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, CA 93940-5121 


Dr. Jovan Lebaric, Code EC/LB....00........c. eee ceecececceeccccceccececeeceees 


Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, CA 93940-5121 


Dr. Clark Robertson, Code ECIRC cecccccccccceccescesesssssessesessteseeuseveseens 


Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, CA 93943-5121 


Kara Kuvvetleri K.ligi MEBS. Bsk.ligi...............essesseeeseees eer 


Bakanliklar, Ankara, Turkey 


Kara Kuvwvetleri K.ligi Teknik ve Proje Yonetim Daire Bsk.ligi ................ 


Bakanliklar, Ankara, Turkey 


Kara Kuvwvetleri K.ligi Kutuphanest......................s::ssesececessseseeeeeceoes 


Bakanliklar, Ankara, Turkey 


Kara Harp Okulu K.ligi Kutuphanesi.......................cccceecsceeeeeeeeeeeees . 


Yenisehir, Ankara, Turkey 


NO YSCl EO Sai asic tis ops eth ees ede else twe sgh aes es decease sapeactennanes: 


Ergenekon Mah. Kumkale Sok. No=27/9 
Yenimahalle, Ankara, Turkey 06210 


105 





No. Copies 


